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Abstract 



Now that is known to be large, a medium baseline reactor experiment can observe the 
fine structure of the electron antineutrino survival probability curve, approximately periodic 
oscillations in L/E with wavelength 47r/|AM|^|. The periodicity with respect to L/E is 
broken by 2-3 oscillations which, in the case of the normal (inverted) hierarchy, shift the 
first 16 oscillations nearly 1% higher (lower) and move the next 16 lower (higher). The 
energy of each peak determines a particular combination of the mass differences, for example 
cos^(6'i2)|AM|]^| +sin^(6'i2)|AM|2| for all peaks visible at baselines under 40 km. Comparing 
these combinations with each other or with NOz/A results one can in principle determine 
the mass hierarchy. Alternately, as the Fourier transforms of the 1-3 and 2-3 oscillation 
probabilities are out of phase by the 1-2 oscillation probability, near the maximum of the 
1-2 oscillation the complex phase of the total survival probability can be used to determine 
the hierarchy. Two interference effects make this task difficult. First, kilometer distances 
between the reactors reduce the amplitudes of peaks below about 4 MeV. Second, even 
reactors 100 or more kilometers away significantly obscure the 1-2 oscillation maximum, 
which also complicates a measurement of the solar mixing angle with a single detector. 
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1 Introduction 



10 years ago Petcov and Piai, on the basis of the existing experimental indications that 
the solar mass splitting AM21 ~ 2.4 x 10~^eV^, suggested that a 20-25 km baseline neutrino 
detector may be able to determine the neutrino mass hierarchy from reactor antineutrinos 
if the solar mass splitting AM|^ is greater than about 10~'*eV^ [1]. It is now known that 
AM21 ~ 7.5 X 10~^eV^ [2], implying that the difference between the normal and inverted 
hierarchies cannot be seen at a baseline below about 40 km. The low fluxes at these long 
baselines led various groups [3, 4] to conclude that the individual peaks in the neutrino 
spectrum would be difficult to resolve, and so the only hope for a hierarchy determination 
would be to sum them via a Fourier transform. Even in this case the necessary detectors 
were extremely large and the experiments slow. 

This all changed with the recent determination of ^13 [5, 6] which, as was hinted last year 
[7, 8], corresponds to a value of sin^ (26*13) up to 10 times higher than had been considered 
by the authors of Refs. [3, 4]. The larger mixing angle increases the sizes of the oscillations 
used to determine the hierarchy by an order of magnitude. This means both that a reactor 
experiment to determine the hierarchy is now practical and also that the analysis of the opti- 
mal baseline and experimental configuration must now be redone. Indeed, such a reanalysis 
is urgent as such an experiment will be built soon [9, 10, 11]. 

We will perform this updated analysis in three papers, mirroring the structure of Ref. [4]. 
In this first paper we will analyze a medium baseline reactor experiment analytically. We 
will derive old observations relating observables of the Fourier transformed spectrum to 
the neutrino mass hierarchy and we will also provide new methods for determining various 
combinations of the neutrino masses and the hierarchy. The optimal analysis will be a 
combination of the old and the new to be determined by simulations. We will also discuss 
problems which have so far escaped attention in the literature. In the second paper [12] we 
discuss the spurious dependence, first observed in Ref. [13], of the parameters of Ref. [4] on 
the theoretical reactor spectrum and on the mass difference | AMI2 1 and we show that it can 
be removed with a weighted Fourier transform. Finally in the third paper [14] will describe 
the results of our simulations, indicating for example the optimal baseline and the optimal 
location for such an experiment in China's Guangdong province. 

A central theme in these papers will be the observation that, since reactors are typically 
separated by of order 1 km in a reactor array, the baselines of neutrinos from different 
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reactors are different. As a result the oscillations of low energy neutrinos, which perform 
half an oscillation while traveling from one reactor to the next, will be invisible at the detector 
as the neutrinos from one reactor arrive at the oscillation maximum while the other is at 
its minimum. We will show that, due to a degeneracy between the hierarchy and a mass 
difference, a determination of the neutrino mass hierarchy is impossible in an experiment 
which cannot resolve these low energy peaks. Thus the angle between the detector and the 
lines extending between reactors is an essential variable in the determination of the optimal 
detector location. For example, for a linear array of reactors like RENO or Daya Bay plus 
Ling Ao, this effect can be eliminated if the detector is placed orthogonal to the array. 

The greatest sensitivity to the hierarchy arises near the 1-2 oscillation maximum, at 
about 60 km. However the flux from a distant reactor at the 1-2 oscillation minimum of 
120 km will dominate over the flux from the near reactor in this region. In fact, this will be 
the case for a detector placed 60 km away from Daya Bay and Ling Ao in the orthogonal 
direction, as the proposed Haifeng and Huizhou reactors will be near the 1-2 minimum. Also 
if a detector is placed equidistant from Daya Bay and Haifeng or Huizhou at the position 
suggested in Ref. [9] or [11] then the 1-2 maximum neutrinos from Daya Bay and Haifeng 
or Huizhou will correspond to the 1-2 minimum for neutrinos from the proposed reactor at 
Lufeng [11]. This not only makes the determination of the hierarchy more difficult, but also 
is detrimental to the sensitivity to 9i2- The ideal solution to this problem is to use two 
detectors at different distances, say 40 and 70 km. However a more economical solution is 
to keep the baselines short so that the flux from the desired reactor complex dominates over 
the fluxes from others, for example one can consider a baseline of 45-50 km. 

We will begin in Sec. 2 with a review of standard results on 3-fiavor neutrino oscillations 
and the electron antineutrino survival probability. We describe the interference between the 
1-3 and 2-3 oscillations [1] which leads to beats at the 1-2 frequency and we numerically 
find the energies of the combined peaks. Then in Sec. 3 we describe how the positions of 
various peaks can be used to obtain combinations of the neutrino mass differences. The 
position of the nth peak determines its own combination AM^^^^. We will see that, as the 
first 10 peaks all determine the same mass combination, they do not allow a determination 
of the mass hierarchy, but the next 5 do, which is why short baselines are not sufficient 
for a determination of the hierarchy. In Sec. 4 we discuss the consequences of the finite 
energy resolution and the finite neutrino flux. Both provide obstacles to locating the n > 10 
peaks, and so to a determination the hierarchy. In Sec. 5 we discuss analyses of the Fourier 
transform of the neutrino spectrum. These have the advantage that, when the nonlinearity 
of the detector response is well understood, they sum the peaks together, and so render the 
signal stronger. We rederive three old ways in which the hierarchy can be determined from 
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this transformed spectrum and add two new methods to the hst. 

Finally in Sec. 6 we discuss the consequences of the fact that not all of the neutrinos 
detected traveled the same distance. The distances may differ by of order a kilometer because 
the individual reactors in an array are not coincident, leading to an interference effect which 
greatly diminishes the amplitudes of the low energy, high n, peaks. We will see that this 
interference can be avoided if the detector is placed perpendicular to the array. Also, while 
it has long been known [3] that neutrinos from reactors at the 1-2 oscillation maximum 
baseline are the most useful for determining the hierarchy, we will see that this signal can be 
overwhelmed by neutrinos from distant reactors at the 1-2 minimum, and we will argue that 
this is indeed the case if a detector is placed at the 1-2 maximum orthogonal to the Daya 
Bay, Ling Ao reactor array. This problem, as well as the related error in a determination of 
6^12, can be reduced by shortening the baseline or, if possible, adding another detector at a 
different baseline. 

Today we received a preprint [13] which has significant overlap with the results in our 
Sec. 2. For example their 0, introduced in Ref. [15], is proportional to our a„/n. 

2 The electron neutrino survival probability 
2.1 Short and long oscillations 

The electron neutrino weak interaction eigenstate is not an energy eigenstate but it 
can be decomposed into a real sum of energy eigenstates 



In the relativistic limit, after traveling a distance L, the survival probability of a coherent 
electron (anti)neutrino wavepacket with energy E can be expressed in terms of the mass 
matrix M 




cos(^i2)cos(^23)|l) +sin(^i2)cos(^i3)|2) + sin(^i3)|3). 



(2.1) 



ee 




(2.2) 
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where AM^^- is the mass squared difference of mass eigenstates i and j. Notice that the 
survival probabihty is a sum of cosines and so its cosine Fourier transform with respect to 
the variable L/E is just a sum of delta functions whereas its sine transform vanishes. 

The three cosines in the survival probability (2.2) identify two characteristic frequencies 
of the L/E oscillations. The P12 term oscillates at a low frequency 

^ ~ 3.8 X 10"^ eV^ ~ 0.19 MeV/km. (2.3) 

Therefore the maximum 1-2 oscillation occurs at 

^ ^ 17km/MeV. (2.4) 



E AMU 2 

A medium baseline reactor experiment, with a baseline of under 100 km, can observe at 
most one or two such oscillations. Instead such experiments will focus on shorter oscillations, 
characterized by the P13 term, which have frequency 



|AM|,| 

r 

2 

corresponding to a wavelength of 



1.2 X 10"^ eV^ ~ 6.1 MeV/km (2.5) 



At a medium baseline reactor one may hope to see 5 to 15 such oscillations. 

What about the P23 term in (2.2)? The frequency is IAMI2I/2, which is about 3% more 
or less than |AM|^|/2 depending on the hierarchy. However the amplitude is less than that 
of Pi3 by a factor of 

a = tan2(^i2) ~ 0.45. (2.7) 

As the frequencies of the two short oscillations are similar but P23 has a smaller amplitude, 
the total short distance oscillation P13 + P23 is a deformation of P13 alone. However the 2-3 
oscillations serve to slightly displace the 1-3 peaks and shift the amplitudes with a pattern 
which repeats at the beat frequency AM|^/2. 

More precisely, while the nth maximum of P13 is at L/E = 47m/|AM|]^|, the nth peak 
of Pi3 + P23 is at 

for a vector a;„ which is determined entirely by the neutrino mass matrix. The positive 
(negative) sign applies to the normal (inverted) hierarchy. As the derivatives of the reactor 
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neutrino spectrum and the 1-2 oscillations are small, the peaks of the total survival proba- 
bility Pee and its product with the reactor neutrino spectrum, which is the observed neutrino 
spectrum, are roughly located at the values given in Eq. (2.8). 

From Eq. (2.8) it is possible to understand the main obstruction to determinations of 
the hierarchy at short baselines. The position of the nth peak is only sensitive to the mass 
combination 

AM^ . = \^MhL_ (2.9) 

As we will see (analytically and then numerically in Fig. 1) at low ra, corresponding to 
peaks visible at relatively short baselines, the values of a„ are roughly linear. Thus in this 
regime the positions of all of the peaks are sensitive to the same effective mass and so are 
independent of the hierarchy so long as that effective mass applies: A change in the hierarchy 
can be compensated by a shift in |AM|]^| while leaving the peaks fixed. The hierarchy can 
only be determined from the nonlinearity of «„. 



2.2 Finding a„ 

The vector encodes the effect of the neutrino mass matrix on the location of the survival 
probability peaks. Therefore a knowledge of this dependence together with a measurement 
of the peaks allows one to reconstruct some elements of the mass matrix. 

The values of a„ are determined from (2.8) by first finding the extrema of P13 + P23 



= A(A3 + P.3)0cA 



cos 1 ' j + tan2(^i2)cos \ 



(2.10) 



cx: sin (±27ra„) + (1 =1= e) tan^(^^i2)sin (27r[±a„ =1= e(n ± a„)]) , e— 



AM|J 



where the upper sign applies to the normal hierarchy. 
As e << 1 and q;„ << we may approximate 

~ sin(27ra„) + tan^(6'i2)sin(27r[a„ - en]). (2.11) 

This can be expanded in a power series in n. The highest energy peaks occur at small values 
of n, where the linear term in this expansion suffices 

~ 27r(l + tan2(^i2))an - 27rtan2(^i2)en (2.12) 

which is easily solved for an 

an ~ esin2(ei2)n ~ O.OlOn. (2.13) 
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Thus the nth peak, for n sufficiently small, lies at 



L Att 



Ann 



(2.14) 



E \AMl 



\AMl\Tsm\9,2)AMl 



where the upper sign corresponds to the normal hierarchy. 

The basic problem facing shorter baseline experiments, which are only sensitive to peaks 
at small n, is already apparent in Eq. (2.14). The mass difference |AM|^| is degenerate with 
the hierarchy 



Therefore any experiment with such a short baseline that all observable peaks have n in 
the regime in which is linear is, alone, incapable of determining the hierarchy. Such an 
experiment can, however, determine the combination 



AM^ff = lAM^il T sm'^{9i2)AM^^ = cos^{9i2)\AMl^\ + sm^{9i2)\AMl^\. (2.16) 



This is just the short baseline effective mass first reported in Ref. [16]. 

Therefore reactor experiments can determine the hierarchy in only two ways. Either 
an accurate determination of AM'^jj can be combined with an accurate determination of 
another combination of the mass differences obtained from another experiment, a difference 
which needs to be known more precisely than the atmospheric mass difference measured by 
MINOS, or else the experiment needs to be sensitive to peaks at large enough n that the 
linear approximation breaks down. So just how large does n need to be? 

2.3 Cubic terms is a„ 

To see where interference between 2-3 and 1-3 oscillations pushes the peaks of P13 + P23 at 
larger values of n, we need to expand to cubic order 



and to substitute this expansion into Eq. (2.11). The linear term in already solves this 
equation at linear order. At cubic order it yields 



AMli I (normal) = jAMlJ (inverted) + 2sin^(6'i2)AM| 



(2.15) 




(2.17) 




(2.18) 



and so 




-1.6 X 10 



-5 



(2.19) 
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n 



0.06 



-0.02 

-0.04^ 

-0.06^ 



The linear approximation to is reliable when the cubic term in Eq. (2.17) is much 
smaller than the linear term 



For example, at the 10th peak the contribution of the cubic term to the energy is only one 
sixth of that of the linear term. The linear term we have seen shifts the effective mass by 
about 1%, thus the cubic term, which is the leading hierarchy- dependent term, only shifts 
the energy of the tenth peak by about one sixth of a percent. The other hierarchy would lead 
to a shift of a sixth percent in the other direction, so overall the difference in the energies of 
the 10th peaks in the normal and inverted hierarchy is about one third of a percent. 

Therefore if only the first ten peaks can be measured at a given baseline, the detector 
will need to be able to determine the position of the tenth peak with a precision of a half 
percent for only a one sigma determination of the hierarchy, making a determination of the 
hierarchy using such an experiment alone quite unlikely. 

The cubic expansion is no longer reliable for higher peaks, but a„ can be determined 
numerically. As seen in Fig. 1, it is periodic modulo 1/e and is zero at every multiple of 
l/2e. But the main problem is that it is nearly linear for n < 10, which is why the hierarchy 
is so nearly degenerate with the shift in the mass differences at all of these peaks. 

K 2/k percent precision measurement of the peak energy of the 14th peak would also 
give a k sigma indication of the hierarchy. The 14th peak can only be seen if it occurs at a 
high enough energy that the neutrino flux and detector resolution are sufficient to discern it. 
For example if one requires it to appear above 3 MeV, so that the detector resolution may 
be better than 2%, then using Eq. (2.6) the minimum baseline is about 



Figure 1: an determined numerically 




(2.20) 



L 



'mm 



45 km. 



(2.21) 
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3 Determining mass differences from peak positions 



3.1 The reactor neutrino energy spectrum 

In the previous section we saw that the aperiodicity of the peaks in the electron (anti)neutrino 
survival probability is determined by the neutrino mass spectrum. In particular, m L/E 
space the peak positions are periodic modulo 1/e ~ 32 and each set of l/2e peaks is displaced 
about 1 percent either left or right depending on the hierarchy. Thus, to determine the 
hierarchy, it suffices to measure the positions of enough peaks to within 2 percent precision. 

The trouble with such a procedure is that no single experiment has access to all of the 
peaks, as there is a single useful basehne per detector and the energy spectrum is limited 
to that which is produced by a nuclear reactor, which leads to a maximum usable neutrino 
energy. Not even all of these energies are accessible as each type of neutrino detector has 
a minimum energy which it is able to detect. To determine which peaks may be seen by a 
particular experiment, one must combine both of these constraints. 

The neutrino flux from a reactor results almost entirely from decays of just 4 isotopes: 
235-[j^ 239p^^ ^^^Pu and ^^'^U. The flux (pii^E) of neutrinos from each isotope i is traditionally 
approximated as the exponential of a polynomial in the neutrino energy E [17] 



Theoretical errors on these fluxes are often claimed to be near the 2-3% level, although recent 
theoretical fluxes [18] appear to be about 6% higher than the fluxes measured at very short 
baseline experiments [19] and at 1 kilometer experiments [20]. 

The precision of such a phenomenological law depends on the degree n of the flt poly- 
nomial. For neutrinos with between 2 and 7.5 MeV, which will be the ones of interest in 
this note, it was shown in Ref. [21] that a quadratic (m = 3) flt tends to introduce errors of 
order 2-3% whereas a 6 parameter (m = 6) flt introduces errors of order 1%, well below the 
error in the theoretical flux. The differences between these parameterizations are oscillations 
over a characteristic scale of 1-2 MeV, which are likely too broad to give a false signal for 
a 1-3 oscillation peak, but may well disguise the depth of 1-2 oscillations and so affect the 
measured value of 9i2. 

The most recent theoretical estimate of the flux is in Ref. [22] , which shows a systematic 
3% excess over the previous estimates [18] at energies above 6 MeV and a 1% surplus beyond 
4 MeV, which is within the theoretical errors of the calculations. This correction again will 
affect a single detector determination of Ou- 




(3.1) 
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E (MeV) 



Figure 2: The theoretical reactor neutrino spectrum as measured with inverse /3 decay. 



Not all of the neutrinos which are generated by a reactor will be measured. The maximum 
number of neutrinos which can be measured at a given energy E is the product of the 
produced flux with the fraction of neutrinos at that energy which can be measured. For 
example, if neutrinos are measured via the inverse /3 decay reaction 



Uf, + p ^ n + e 



(3.2) 



then the maximum number of neutrinos detected is the flux/area at the baseline L multiplied 
by the inverse /3 decay cross section which at tree level is [23] 



a{E) = 0.0952 x lO^^^cm^E^pjUeV 



(3.3) 



where the positron energy and momenta are, ignoring the neutron recoil, given in terms of 
the neutron, proton and electron rest masses 



Ee = E -mn + mp + mer^ E - 780 keV, pe = ^ El - 



mi 



(3.4) 



Combining the reactor flux (3.1), with the coefficients a^j and typical isotope fractions 
of Ref. [17], with the tree level cross section (3.3), one obtains the theoretical reactor flux 
depicted in Fig. 2. While inverse /3 decay is kinematically forbidden \i E < m„ + me — vHp ~ 
1.8 MeV, it can be seen in Fig. 2 that the fiux/energy is maximized at 3.6 MeV, falls to 
one third of its maximum by 6 MeV and to 10% of its maximum by 7.5 MeV. Thus useful 
information can only be obtained about the spectrum for energies within a factor of 2, which 
for each detector corresponds to values oi L/E and so n within about a factor of 2. 
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Figure 3: Theoretical neutrino fluxes, including 3 flavor oscillation, for the normal (solid 
curve) and inverted (dashed curve) hierarchies as seen at 40, 50, 60 and 70 km. 

3.2 Effective masses at various baselines 

When neutrino oscillations are included, the theoretical flux $(£') from a reactor at a distance 
L is then multiplied by the survival probability Pee given in Eq. (2.2) 

$(E) = cME)cr{E)PUL/E) (3.5) 

i 

where is the quantity of each isotope in the reactor. 

We will use AM|^ and sin^(2^i2) from Ref. [2], IAMI2I determined by combining neutrino 
and antineutrino mass differences from Ref. [24] and sin^(26'i3) from [25] 

AMa^ = 7.50x10"^ eV^ lAM^al = 2.41x10'^ eV^ sin2(2^i2) = 0.857, sin2(2ei3) = 0.089. 

(3.6) 

|AM|]^| is determined using the normal and inverted hierarchies. The resulting neutrino flux 
is shown in Fig. 3 at baselines of 40, 50, 60 and 70 km. 

The numbers n of the local maxima can be read from (2.8) by approximating |AM|^| ~ 
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Figure 4: Theoretical neutrino fluxes, including 3 flavor oscillation, at 40 and 58 km for both 
hierarchies with the same value of AM^jf = 2.46 x lO'^ eV^ In the left panel, at 40 km, 
the hierarchies are difficult to distinguish because a„ is nearly linear at the visible peaks. 
This degeneracy is broken by the higher n peaks visible at 58 km, seen in the right panel. 



IAM32I and setting = 

n ^ — ~ 0.97-— — — — . 3.7 

Att E E/MeV ^ ' 

As the error on IAMI2I is of order 3% and the difference between AM^j-j and jAMfgl is also 
of order 2%, one may expect an error of 3-5% in Eq. (3.7). The fractional energy difference 
between the nth and {n + l)st peak is 1/n. Therefore an optimal detector can robustly 
determine n given a single peak only if n is less than about 10. Limitations placed by the 
finite energy resolution of the detector are discussed in Subsec. 4.1. 

Although the energies of peaks in Fig. 3 are strongly dependent upon the hierarchy at 
every baseline, this does not mean that a measurement of the spectra can actually allow one 
to determine the hierarchy. The problem, as was discovered in Ref. [15] in this context and 
as we will describe in Sec. 2, is that the energies of the first 10 or so peaks only determine 
the mass difference AM^^jr of Eq (2.16), leaving the hierarchy degenerate with |AM|^|. 

Thus if the baseline is short enough so that only these peaks may be reliably measured, 
then there will be a value of AM^^j that reproduces the peaks for both hierarchies, as seen 
at 40 km in the first panel of Fig. 4. Here the peak n = Q can barely be seen at 6.6 MeV, 
whereas n = 7 at 5.7 MeV is clearly discernible. The largest peaks are n = 8, 9, 10, and 11. 
However the energies of these peaks are independent of the hierarchy at constant AM^j^. 
The hierarchy-dependence becomes somewhat larger at the 12th peak, which is located at 
about 3.32 MeV in the case of the normal hierarchy and 3.30 MeV in the case of the inverted 
hierarchy. Thus if AM^^j peaks can be determined at better than 1% from the low n peaks 
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and then the location of the 12th peak can be determined with a precision of better than 
1%, a determination of the hierarchy would be barely possible. The lower energy peaks are 
smaller, but more hierarchy dependent. For example, the 16th peak would be at 2.50 MeV 
with the normal hierarchy, but 2.46 MeV with the inverted hierarchy. This 2% difference is 
just within the resolution of the proposed detectors of Refs. [4, 9, 11], and so when combined 
with an accurate measurement of AM^j^j one could potentially determine the hierarchy at 
the 1-2(7 level. 

In the second panel one can see the electron neutrino survival probability at 58 km with 
both hierarchies and the same value of AM^jjr. At this long baseline one can see maxima up 
to n = 20, where the nonlinearity in a„ is appreciable and so, as one can see, the low and 
mid energy peaks are hierarchy-dependent at the 1-2% level. 



3.3 The 1-2 oscillation maximum 

It is clear from Fig. 4 that if AM^jy is determined from the low n peaks then the locations 
of the peaks at the 1-2 oscillation maximum 

n = — ~ 16, E= — ^MeV (3.8) 
2e 18 km ^ ^ 

are hierarchy-dependent, and so one may hope to use their locations to determine the hier- 
archy. This can be done if the peak locations allow for a combination of the mass differences 
which is distinct from AM^^j. In fact, two such combinations can be determined, one from 
the location of the peaks and one from the distance between them. 

To derive these two combinations, we will need to find a„ near the 1-2 oscillation maxi- 
mum. This is easily obtained by expanding (2.11) about n = l/2e 



~ sin(27ra„) — tan^(6'i2)sin ( 27r 
Linearly expanding the sine function yields 



1 

a„ — e \ n 

2e 



(3.9) 



0~(l-tan^(^i2))an + e(^-^ ) (3.10) 



and so for n ~ l/2e 



~ — ^^hr^- (3-11) 



tan^i^(7i2 



At n = l/2e, a„ = 0. Of course n must be an integer, so it can never be precisely equal 
to l/2e. Nonetheless, a„ will be of order 0.01 when n is at the closest integral value, leading 
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to a less than 0.1% contribution to the energy. When = 0, the energy of the peak is 

inriL 2t[L 
" = WmM = AMI 

allowing for a precise determination of AM|^. Of course, to know that one is at the 1-2 
maximum by counting peaks, one must know e and so this is related to a measurement of 
|AM|]^|. Alternately, one can find the 1-2 maximum by looking at the large oscillations in 
the flux due to 1-2 mixing and then count peaks to determine n at the minimum, which 
allows for a determination of e and so |AM|^| from AM|]^. 

The distance between the peaks near the 1-2 maximum allows for a determination of the 
mass differences which is independent of precise knowledge of the 1-2 oscillation parameters. 
Inserting the linear expansion of a„ (3.11) into the definition of a„ (2.8) one finds that the 
distance between two peaks is 

L\ 47r 

fail V i - tan-|^t/i2 

and so it determines the effective mass 



AML. = ( 1 ± ^..T^ 1 lAM^ I = 2 t^^'(^i2) |^^2 I L^iAM^ 



V l-tan2(ei2)y l-tan2(^i2) l-tan\9i2)^ 

(3.14) 

This effective mass is quite different from that defined in Eq. (2.16). Approximating 
tan^(^i2) = 1/2 one finds that the spacing between the first 10 peaks yields an effective mass 



AM^ff = l\AMl\ + l\AMi,\ (3.15) 



while the peak spacing between n = 14 and n = 18 yields 

AMl^ = 3\AMl\ - 2\AMi,\. (3.16) 

A detector at a baseline of less than 50 km can accurately determine the combination (3.15) 
while one with a baseline between 45 and 70 km, as it sees higher n peaks, may more easily 
measure the combination (3.16). The normal mass hierarchy is equivalent to the second 
mass being larger than the first, and so by comparing these masses one can determine the 
hierarchy. 

The difference between these masses is quite large, about 7%, however a 7% measurement 
of AM^j„ requires a 7% precision in the measurement of the difference between the the peaks 
in the linear regime near to the 1-2 oscillation maximum n = l/2e. As can be seen in Fig 
1, this linear regime includes about 8 peaks, corresponding to a 40% variation in energy. 
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This means that a 7% precision in the distance between the peaks requires a 4% precision 
in the energy of the peaks, much less than was required at shorter basehnes. Therefore a 
comparison of the low n peak positions and the 1-2 maximum peak separations is a promising 
test of the hierarchy, so long as the statistics are sufficient for the peaks to be observed. 

4 Resolution and flux constraints 
4.1 Energy resolution 

The neutrino energy E is determined from the positron energy by adding 780 keV (3.4). 
The positron energy is determined by counting photoelectrons in a scintillator. The number 
of photoelectrons is proportional to the positron energy, and so the energy resolution cr^; is 
proportional to the square root of the positron energy. For example, at Daya Bay II it has 
been suggested in Ref. [4] that a resolution of 



sume this resolution in what follows. The observed positron energy spectrum is then the 
convolution of the true spectrum with a Gaussian of width a 



Both spectra are displayed in Fig. 5. 

The observed spectrum Pobs(-£'e) is the best that can be hoped for, once the response 
of the detector is understood and with an infinite number of neutrinos. Finite fiux effects 
will be briefly discussed in Subsec. 4.2 and then discussed in detail in the companion paper 
on our simulation results. A poorly understood nonlinear detector response is fatal to the 
Fourier analysis that will be discussed in Sec. 5, but individual peaks may be analyzed where 
the response is understood. 

Even in this idealized setting, it is clear from Fig. 5 that the low energy (high n) peaks 
cannot be resolved. How high is n for the lowest energy peak that can be resolved? This 
depends on the baseline and the neutrino flux. However a rough answer is obtained by 
asserting that the distance between the nth maximum and the adjacent minimum (2.6) 



aE = 0.03^/{Ee)MeV 



(4.1) 



will be achieved, corresponding to a fractional resolution of 3% / E /MeV . We will as- 
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Figure 5: The true neutrino spectrum P and the observed spectrum Pobs with a resolution 
of 3%/a/E' using the normal hierarchy at baselines of 40, 50, 60 and 70 km. Notice that 
the lowest energy oscillations are smeared away and the energy threshold for this smearing 
increases with the baseline, such that at higher baselines the maximum n observable increases 
slowly. 
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be greater than 



2a = 0.06x/(^e)MeV = 0.06 MeV^^O-g?^^^^ - 0.78. 



(4.4) 



These two quantities are equal when 



L/km = 7.4 X 10" V ± ^ 1 - ^ j . (4.5) 

The equahty with the plus sign yields the upper bound on observable n for a given baseline 
L, whereas the expression with the minus sign typically yields a value of n at energies where 
no reactor neutrinos are observed. 

Alternately (4.5) provides the baseline L necessary to observe the nth peak. When n < 15 
it provides no bound at all, so long as there is enough flux and the detector response is well 
understood, the energy resolution is not an obstruction to observing these peaks. Of course 
at low baselines they may be difficult to observe because there simply are not many or any 
neutrinos observed at the corresponding energy. At n = 16, 17 and 18 one finds minimum 
baselines of 45 km, 57 km and 70 km respectively. In practice the true minimum depends 
on the neutrino fiux. But this rough estimate shows an essential point, that with a 3%/\/E 
fractional resolution any medium baseline, between about 40 and 70 km, is sufficient to 
observe the 1-2 oscillation maximum if there is enough neutrino flux. Longer baselines only 
marginally extend the reach to higher peaks, although each of these peaks is in the 1-2 
maximum region and so even just 1 or 2 more peaks can greatly enhance the possibility of 
determining the mass hierarchy. 

While the maximum observable n is reasonably independent of the baseline, this deriva- 
tion shows that it is strongly dependent upon aE- A resolution comparable to that of Daya 
Bay or RENO would lead to a maximum n which is still in the linear region of a„, and so 
the hierarchy would be unobservable. 



4.2 How much neutrino flux is required? 

The neutrino flux that can be observed by a large, distant detector via inverse /5 decay is still 
quite unknown, even the efficiency of the detector is difficult to predict. A rough estimate 
can be made using the flux observed at Daya Bay, using the flux normalization-independent 
determination of ^13 to eliminate the loss due to 1-3 operation. At Daya Bay 17.4 GW of 
thermal power yields 80 neutrinos/day at each 20 ton detector at a weighted baseline of 1600 
meters. Therefore the total measured flux/year from a P GW reactor complex measured at 



16 



a detector of mass M at a baseline L is roughly 

P M 2.6 km^ ,(P/GW)(M/ktons) , , 

$ = 365 x 80 x— — — -— = 2.2 x 10^ ^ ' , ' ^. 4.6 

17.4 GW 0.02 ktons (L/km)^ ^ ^ 

For example, at a distance L from the Daya Bay and Ling Ao reactors a 20 kton detector 
may observe 7.6 x 10''/L^ neutrinos/year, where the length L is measured in kilometers. 

Let the energy width of a peak be AE, the integrated flux within that energy range be 
(f) and the peak be a fraction A higher than the flux in that energy range from the reactor in 
question considering 1-2 oscillations and also the flux from distant reactors. In the absence 
of other reactors the fraction A varies from sin (26*13) at low n to 1.7sin2(26'i3) at the 1-2 
maximum. The fractional error in the flux in that range will be 1/ a/0. Therefore the peak 
may be observed if ^/(j)A > 1. 

Simply observing the peaks is useful for two reasons. First of all, the distance between 
peaks is roughly L|AM|^|/7r and so by identifying peaks, one can check the consistency of 
the location of other peaks. Second, recall that it is easier to determine the n value of the 
well-separated, high energy, low n peaks. If one can observe the peaks from low to high n 
then it is possible to count them and so determine the n values of the low energy peaks. 
While the hierarchy can be determined from the distance between the high n peaks, as 
described in the previous subsection, without knowing the precise value of n, nonetheless if 
one knows the n value of a high n peak it can be used to determine a combination of the 
mass differences via (2.8). This determination has a precision of 1/n, and so at high n it 
leads to a precise determination. 

However, to determine the hierarchy it is not enough to observe the peaks, one must 
determine their energies as precisely as possible. How precisely may they be determined? 
They may be determined within an energy 6E if the neutrino surplus can be seen in width 
6E bands within the peak. This requires 

1 < = V^Afn (4.7) 

where / is the fractional energy precision desired. This implies that the maximal fractional 
precision with which the energy of a peak can be measured by a detector with perfect 
resolution is 



As mentioned above, with no unwanted backgrounds from other reactors, A varies between 
0.09 at the 1-2 maxima to 0.16 at the minima. 

Consider for example the 11th peak at a 40 km baseline, which lies at 3.6 MeV. At this 
point A ~ 0.12. The flux within the peak is about 4% of the total flux (4.6), which each 
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year at a 20 kton detector at 40 km from Daya Bay may be 



(j) = 0.04 X 7.6 X 107(40)2 = 1.9 x 10\ (4.9) 
Therefore, after m years, the best resolution of an ideal detector would be 

■^'""^ " 0.12 X lOv^m X 103 " 577^' ^^'^^^ 

Thus an ideal detector can find the 11th peak energy to within l/i/m percent after 2.5m 
years. As was described in Subsec. 4.1, since the peak width is greater than the high energy 
resolution proposed at a new medium baseline detector experiment in Ref. [4], resolution 
effects will not significantly alter the determination of this peak. Therefore one may expect 
AMgj.j to be determined to a precision better than 1 percent at a 40 kilometer baseline 
experiment with no backgrounds from other reactors. 

However, to determine the hierarchy, one also needs to determine another combination of 
the neutrino mass differences. This requires the measurement of a higher n peak. Consider 
for example the n = 16 peak at 2.5 MeV. While it is in general a poor approximation to 
ignore the backgrounds provided by distant reactors, if one does ignore them then, since the 
peak is near the 1-2 maximum, the relative peak height is A = 0.16. The total flux in the 
peak is only about 0.8% of the total, and so 

1 1 

" 0.16 X 16^3:8771 X 102 " 5V^' ^^'^^^ 

As can be seen in Fig 4, a 3% precision is sufficient to determine the hierarchy and so the 
energy of this peak would yield a 2a determination after 2 to 3 years. As the width of the 
peak is about 3%, a 3%/a/E resolution may reduce the amplitude by a factor of 2, still 
allowing for a determination of the hierarchy. 



5 The Fourier transform of the survival probabiUty 

While each peak provides some information regarding a combination of neutrino mass 
differences and therefore the hierarchy, it may well be that the fluxes are too weak or the 
backgrounds too large for the peaks to be reasonably well identified. Complimentary infor- 
mation can be obtained by combining the peaks. As the electron survival probability Pee is a 
sum of periodic cosine functions, they can be combined by a Fourier transform. Even when 
individual peaks are hard to identify, the combination probed by the Fourier transform may 
well be visible and so may provide the best chance for determining the hierarchy [3] . 
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5.1 The complex Fourier transform 



For simplicity we will approximate the observed electron antineutrino spectrum by a Gaus- 
sian distribution in L/E space 

=e-(^-^<4))V (5.1) 

where (l/E) is the average 1/E of a detected neutrino. While with only slightly more 
complicated equations the following could be avoided, we will make the crude approximation 
that (5.1) is the neutrino spectrum after 1-2 neutrino oscillations, and that the expectation 
value of the inverse energy is therefore taken with respect to the 1-2 oscillated spectrum, 
which depends upon L. 

1-3 oscillations affect this spectrum by introducing a modulation equal to (^{L/ E)Pi^{L/ E). 
The Fourier transform of this modulation is 

Fi3(fc) = I d(^\ ^{L/E)P,,{L/E)e^'f (5. 



where we have factored the sin^(2^i3) out of the definition of Ei^. This quantity has two 
peaks, one ai k = \AM^^\/2 and one at k = —\AM^-^\/2. At each peak, one of the two terms 
in the parenthesis dominates, and the other is suppressed by a factor of order e"^/^, which 
is large enough that we will neglect the sub dominant term. Therefore, near the first peak 

^ ^^W(M,-* ('-^)',<-^)^<i'. (5.3) 

The norm of the complex function F13 has a maximum at = |AM|j^|/2, where F13 is 
real. The real part of -F13, corresponding to a cosine transform, is, within the validity of 
the approximations described above, symmetric about this maximum. The imaginary part, 
corresponding to a sine transform, vanishes at this maximum and is antisymmetric about it. 

The transform E23{k) can be calculated identically. Near the positive k maximum the 
sum of the two transforms is just 



a\/Tc 



crwn 
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where the + sign apphes to the normal hierarchy. 

The first term corresponds to the old peak, at k = |AM|^|/2 and the second to a new 
peak at k = IAMI2I/2, which on its own would have a height equal to tan^(^i2) ~ 1/2 times 
that of the first. Of course, due to the phase difference between the two terms, for some 
choice of parameters the interference between the two terms implies that the second is not 
a local maximum of the norm of the Fourier transform. The k value of the peak then gives 
the corresponding mass difference. Thus if the smaller peak is to the left, corresponding to 
smaller k, of the larger peak then IAMI2I < |AM|^| and so one can conclude that there is a 
normal neutrino mass hierarchy [3]. 



5.2 Determining IAM3J at the 1-2 oscillation maximum 



We will refer to the baseline 



271 



(5.5) 



as the 1-2 oscillation maximum, as it is roughly the baseline at which the largest fraction 
of neutrinos disappears as a result of P12. It is approximately 58 km. At this distance the 
Fourier transform (5.4) simplifies as the two terms are precisely out of phase 



awn 



cos2(^i2)e ^ 



3ll\ ( 

- sin^(6'i2)e 



|AJW|2|^ ^' 



(5.6) 



The cosine transform is just the real part of F. Up to /c-independent factors, it is 
proportional to 



FU~k)= {e-'r>'-tan\e,,)e~' 



COS 



2'Kk 
AMI 



(5.7) 



where we have defined the distance from the 1 — 3 peak to be 



(5. 



The maxima of Fcos are found by setting to zero its derivative with respect to k, yielding 
the condition 
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where again the + sign corresponds to the normal hierarchy and the — sign to the inverted 
hierarchy. The largest peak is the closest to = 0, the absolute maximum of the Fourier 
transform of P13. To find this peak we may expand k about 0, at linear order we find 

1 + g {e^cot\eu) - 1) + 2/3 



where we have defined the constant 



(5.11) 

lb 



As the denominator of (5.10) is much greater than 1, we learn that 



\k\ « (5.12) 



and so the maximum of Fcos lies at approximately 

I AM. 



This means that if the detector is placed at the 1-2 oscillation maximum baseline then the 
peak of the cosine transform of the full neutrino spectrum lies essentially at the peak of P13 
alone. This is not because because the frequencies of the P13 and P23 terms are similar, but 
because as seen in Fig. 6, at this baseline the absolute maximum of the Fourier transform of 
Pi3, which corresponds to its frequency, happens to coincide with one of the minima of the 
cosine transform of P23, which is not its frequency. Thus, at the 1-2 oscillation maxi- 
mum, the absolute maximum of the sum of two cosine transforms is coincident 
with that of P13 alone, allowing for a direct and precise determination of |AM|]^|. 
This may be useful on its own even if the hierarchy has already been determined by NOvA 
or T2K. 



5.3 Determining the hierarchy at the 1-2 oscillation maximum 

Can this simplification also yield information about the hierarchy? A precise determination 
of |AM|J combined with MINOS' best fit for lAMfal could give a la answer to this question, 
within 5 years MINOS+ may improve this to 2a. But with enough flux the peak structure 
of the Fourier transform alone yields some information about the hierarchy. 

As can be seen in Fig. 6, just as the relative sign in (5.6) implied that the global maximum 
of the cosine transform of the total spectrum is essentially coincident with that of P13, it 
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Figure 6: The cosine Fourier transform of P12 (blue), P23 (purple), P13 (green) and Pee 
(yellow) at 58 km in the case of the normal hierarchy. Note that the P13 and P23 curves are 
one half of a wavelength out of phase, so that the total extrema coincide with those of P13. 



also implies that the minima just to its left and right are roughly coincident with those of 
L _ |AM|i|-AM|i J, _ \AMl\+AMl 

"'mill 2 ' min 2 ' W'-'-^J 

These are minima for the simple reason that the cosine on the right of (5.7) is equal to —1. 
Substituting these values of k into (5.7) we find that the values of the cosine transforms at 
the two minima are 

(5.15) 

where as usual the upper sign corresponds to the normal hierarchy. In the case of the normal 
hierarchy, the positive sign means that the second term is larger on the left, which implies 
that the minimum on the right is deeper. In the case of the inverted hierarchy these two 
depths are interchanged, and so the peak on the left is deeper. Thus we have recovered the 
criterion for determining the hierarchy from the cosine transform which was proposed in Ref. 
[4]. 

A similar analysis can be applied to the imaginary part of F{k), which is obtained via 
a sine transform. As P(|AM|]^|/2) is real, the sine transform vanishes at the maximum of 
the global cosine transform. The sine transform then has a maximum on the right and a 
minimum on the left. The opposition of the phases of the Fourier transforms of P13 and P23 
at the 1-2 oscillation maxima again imply that these extrema of the sine transform of the 
full spectrum roughly coincide with the extrema of the sine transform of P13 alone, as can 
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Figure 7: The sine Fourier transform of P12 (blue), P23 (purple), P13 (green) and Pee (yellow) 
at 58 km in the case of the normal hierarchy. Note that the P13 and P23 curves are one half 
of a wavelength out of phase, so that the total extrema coincide with those of P13. 



be seen at 58 km in Fig. 7. They simply correspond to the values of k for which the phase 
in (5.6) is ±n/2 
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Defining the sine transform so as to take the imaginary part of F with the same normalization 
as for the cosine transform one then finds 
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(5.17) 

where the upper sign corresponds to the normal hierarchy. In the case of the normal hierarchy 
the second term in the minimum on the left is larger and so |Psin(^^)| < |-^sin(^'^)|5 whereas 
in the case of the inverted hierarchy the maximum on the right is larger than the minimum 
on the left. Thus we have reproduced the correlation between the hierarchy and the relative 
sizes of these extrema observed in Ref. [4]. 

Note that in the case of the normal hierarchy both extrema are more positive and in the 
case of the inverted hierarchy both are more negative, in other words near its maximum 
F has a positive (negative) imaginary part in the case of a normal (inverted) 
hierarchy. This provides a new test which allows one to extrapolate the hierarchy from 
the Fourier transform of the survival probability. Of course the optimal indicator of the 
hierarchy will be a weighted sum of these three tests, with weights that may be determined 
by series of simulations. The ability to distinguish the hierarchies may also be improved by 



23 



convoluting the observed spectrum with a slowly varying function of {L/E) which weights 
neutrinos near the 1-2 maximum more heavily. One can also use simulated data to determine 
which weighting functions are optimal for this task. 



5.4 Nonlinear Fourier transform 

The Fourier transform methods described above essentially work because the phase at the 
maximum is determined by the hierarchy, positive for the normal hierarchy and negative for 
the inverted hierarchy. In the case of just 1 — 3 oscillations, this phase is zero, since the cor- 
responding oscillations are a cosine and the real part of the Fourier transform is determined 
by the cosine transform. However when P23 is included the peaks of the untransformed spec- 
trum move according to Eq. (2.8). The distance between the untransformed peaks changes, 
which in general moves the Fourier transformed peak. 

Critically, the untransformed P13+P23 also loses its mod 2/\AM^^ \ periodicity, as a„ is not 
a linear function of n. A given detector is only sensitive to some of the peaks, corresponding 
to a certain region in Fig. 1. Such regions are generally dominated by a domain in which 
an can be approximated not by a linear function, but by a linear function plus a constant 
offset. This constant offset implies that the convolution of P13 + P23 with the reactor neutrino 
spectrum is not of the form cos{k,L / E) but rather of the form cos{k,L/E + c) where the sign 
of c is determined by the hierarchy. This offset, c, leads to a translation in L/E space which, 
after the Fourier transform, becomes a phase. Thus the hierarchy determines an overall 
phase of the Fourier transform, which leads to the observable indicators described in the 
previous subsection. 

The Fourier transform method is robust since it sums multiple peaks together, and so 
it requires less neutrinos than a direct analysis of the positions of the peaks. However it is 
inefficient because, as was just described, it works by approximating the function a„ to be 
affine in the energy range which is probed. So one might ask if the performance would be 
increased by performing not an ordinary Fourier transform, but a Fourier transform with 
the nonlinearity of a„ built in. The nonlinearity depends upon n ± a„ which depends on 
the hierarchy and also weakly upon the neutrino mass matrix. One may therefore attempt a 
nonlinear Fourier transform with both choices of hierarchy, and if desired, a weight function 
g{L/E). Such a nonlinear cosine transform is given by 

m = I d(^^^^L/E)PUL/E)g{L/E)cos(^k^T2na(^^^^^ (5.18) 

where a{n) is a possibly hierarchy-dependent function which interpolates between the dis- 
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Crete values of and the positive sign corresponds to the inverted hierarchy, such as 
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;i ± e)cot2(6'i2) + cos(27ren) 



(5.19) 



This transform will add all of the peaks together with the same phase for the correct hierarchy 
whereas the peaks will be distorted by the other hierarchy. Therefore the correct hierarchy 
can be determined from the fact that the corresponding nonlinear Fourier transform will 
have a larger absolute maximum. 



6 Interference effects 

6.1 Interference between reactors separated by 1 km 

Consider two reactors of equal thermal power separated by a distance D. If a detector is 
located at a distance L » D from the nearest and if the baseline makes an angle 6 with 
respect to the line passing through both reactors, then the distance from the detector to the 
far reactor will he L + d where d = Dcos{6). 

Adding the flux from both reactors, one finds that the 1-3 oscillations interfere 

P.3 0. COS f^^^V- (^^^^) = 2C0S 0-^) COS n-m\(^y 
\ 2E J \ 2E J \ AE J \ 2E 

(6.1) 

The neutrinos from the two sources are not coherent, it is not the wavefunctions that add, 
but the probabilities. And the result is that the amplitude of the oscillations at energy E is 
suppressed by a factor of 



\AMl\d\ ( rf/km \ 



In particular they annihilate entirely when 

A = 0.5-^. (6.3) 
km MeV ^ ' 

Consider for example the pair of Daya Bay reactors and the two pairs of Ling Ao reactors 
which lie along a line at a distance of 0.8 and 1.4 km from the Daya Bay reactors. The Daya 
Bay II reactor location suggested in Ref. [9] is at an angle of 15 degrees with respect to 
a continuation of this line, yielding d =1.2 km for the nearest and furthest reactor pair. 
The resulting change in the spectrum is shown in Fig. 8. The 1-3 oscillations at 2.4 MeV 
completely cancel between these two reactor pairs, leaving only the oscillations of the middle 
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Figure 8: The spectrum observed at the Daya Bay location suggested in Ref. [9] assuming 
the normal hierarchy. The solid black curve includes the interference between Daya Bay and 
Ling Ao reactors and the red dashed curve does not. It is assumed that the layout of the 
Haifeng reactor complex yields as much interference as that at Daya Bay and Ling Ao. 

reactor, and so effectively damping the oscillation amplitude by a factor of 3, which implies 
that an equally precise measurement of a peak at that energy requires 9 times as much flux 
as it would have without interference. 

At the peak energy of 3.6 MeV the annihilation is not complete, but the 1-3 oscillation 
amplitude of neutrinos coming from the nearest and farthest detectors is reduced by a factor 
of cos(l) = 0.6. Thus the total amplitude of the peak from all three reactor pairs is reduced 
by about 30%, and so 70% more neutrino flux will be necessary to observe these peaks. 

This problem can be avoided if the detector is equidistant from all of the reactors. In the 
case of the Daya Bay and Ling Ao complex, in the case of RENO and somewhat trivially 
in the case of Double Chooz this is possible as all of the reactors essentially lie along a line. 
It means however that such a detector will not be equidistant from any reactor that may 
eventually be built at Haifeng or Huizhou, thus reducing the neutrino fluxes assumed in 
Refs. [4, 9, 11] by a factor of 2. Worse yet, the Haifeng and Huizhou reactors then provide 
a strong and undesirable background, as will be described in Subsec. 6.2. 
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6.2 Interference between reactors separated by 100 km 



The 1-2 oscillation maximum 
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(6.4) 



provides an ideal baseline to determine the mass hierarchy for a number of reasons. Among 
these is that while the P13 + P23 oscillation amplitude of (cos^(^^i2) — sin^(^^i2))sin^ (2^^1,3) is a 
factor of 3 less than its amplitude sin^(2^i3) at the 1-2 maxima, the flux remaining after 1-2 
oscillations is also smaller by a factor of 1/(1 — sin^(^i2)cos^(6'i3)) which is about 5. Thus 
the 1-3 oscillations are a larger fraction of the total flux, making them easier to see above 
systematic errors, although not above statistical errors as 3 > -\/5. 

However the smaller signal and smaller flux means that this part of the spectrum is 
particularly prone to interference from distant reactors, in particular those that are near the 
flrst 1-2 oscillation minimum 



In this case the additional factor of 5 in flux from the distant reactor outweighs the factor 
of 4 distance suppression, and so if both reactor complexes have the same strength then 
most of the flux at the 1-2 maximum will be background. This means that to obtain the 
same energy resolution at the 1-2 maximum peaks one will need more than 4 times as much 
neutrino flux. 

This is the case for example with a detector placed 58 km away from the Daya Bay /Ling 
Ao complex, perpendicular to the reactors. It would be at the 1-2 minimum of the proposed 
Haifeng and Huizhou reactors and would also suffer signiflcant contamination from the Tais- 
han and Yangjiang reactor complexes, at each of which at least 3 reactors are already under 
construction. Similarly a reactor placed 60 km from Daya Bay, Ling Ao and Haifeng would 
be at the 1-2 minimum of the proposed 17.4 GW thermal power reactor complex at Lufeng. 

As this effect increases the flux at the 1-2 oscillation maximum, it is also a serious obstacle 
to an accurate measurement of 9i2- If a model of the background neutrino flux from distant 
reactors is wrong, the error in the total flux can be compensated for by an error in the 
determination of 9i2- Ideally this problem can be solved by using two medium baseline 
detectors instead of one, which would break this degeneracy. However in a world limited 
by costs, it may instead be necessary to simply try to make the background from other 
reactors as small as possible, thus minimizing the potential error in the determination of the 
hierarchy and in the determination of 9i2- 




(6.5) 
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Unlike the short distance interference problem discussed in Subsec. 6.1, the fractional 
contamination caused by distant reactors depends on the baseline L to the reactor whose 
neutrinos provide the signal. The fractional contamination from undesired reactor neutrinos 
is inversely proportional to the desired signal strength, and so it is proportional to . Thus 
this problem is minimized by placing the detector as near to the reactor as possible. If there 
is only one detector, and one wishes to use it to determine the hierarchy, then as explained 
in Subsec. 2.3, this minimum distance cannot be shorter than about 45 km. 

7 Conclusions 

The newly discovered high value of ^13 means that the determination of the neutrino 
mass hierarchy at a medium baseline reactor experiment is now practical. Previous analyses 
of such experiments have assumed values of sin^(2^i3) an order of magnitude or more below 
its true value. At such low values, individual peaks of the electron antineutrino spectra could 
not be resolved and one instead needed to rely upon a Fourier analysis [3, 4]. 

In this note we have reconsidered the determination of the neutrino mass hierarchy now 
that 1-3 oscillations are large and so individual 1-3 peaks in the spectrum can easily be 
observed. We found that the position of each peak determines a particular combination 
of the neutrino mass differences, for example the first 10 peaks all determine the same 
difference cos {O12 sin (^12) 1^-^32 1- particular this means that the first 10 peaks 

alone cannot be used to determine the hierarchy, as they are fit equally well by the wrong 
hierarchy model in which this mass difference is preserved. On the other hand we found that 
the positions of the next 10 peaks depend on distinct combinations of the mass differences, 
and so combining the energies at different peaks with some beyond the 10th can lead to a 
determination of the hierarchy. We also estimated the detector resolution and neutrino flux 
which are needed for such a goal, although an accurate determination will be left for the 
simulations to be discussed in our companion paper. 

The information about the hierarchy is therefore contained in the low energy peaks, 
which suffer the most from the effects of poor resolution, low neutrino flux and interference. 
While the individual peaks are somewhat difficult to resolve, the situation can nonetheless 
be improved with a Fourier transform. We have analyzed the Fourier transform of the 
spectrum, deriving the phenomenological hierarchy indicators proposed in Refs. [3, 4] and 
providing a new indicator, the complex phase at the peak of the Fourier transform, which 
we claim will be positive for the normal hierarchy and negative for the inverse hierarchy. 
We also propose a new hierarchy-dependent nonlinear Fourier transform which will lead to 
a higher peak using the transform corresponding to the correct hierarchy. Our simulations 
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show that an analysis which combines the hnear and nonhnear transform methods is able 
to identify the hierarchy considerably more reliably than one which relies entirely upon the 
linear transform. 

The strength of the Fourier transform method is that it sums together all of the peaks 
to increase the strength of the signal with respect to the noise. If the energy spectrum 
is shifted, this simply leads to an overall phase in the Fourier transform which does not 
seriously affect the analysis. Likewise a uniform stretching of the spectrum simply shifts 
the peaks, which does not affect the determination of the hierarchy at all. Any nonlinearity 
however poses a much more serious problem, as it can lead to an interference between the 
peaks in the spectrum which mutates or destroys the peak structure of the Fourier transform. 
If the nature of the nonlinearity is known then one can adjust the analysis to correct for 
it [26] however if the nonlinearity is known it could be corrected directly in the reading 
of the energy. In general the nonlinearity of the response is only known at energies where 
the detector has been calibrated with radioactive sources. It may therefore be desirable to 
modify the Fourier transform so as to weigh the more reliable energies more heavily. One 
may also wish to weigh the energies near the 1-2 maximum more heavily, as they contribute 
few neutrinos but are indispensable in a determination of the hierarchy. These weights can 
be optimized by testing various hierarchy determination algorithms against simulated data. 

We also discussed two interference effects. First of all, reactors in the same array are 
generally separated by distances of order 1 km, which means that neutrinos arriving at a 
detector from one reactor at a 1-3 maximum may be at the same energy as those arriving 
from another at a 1-3 minimum. The result is that the 1-3 oscillation signal can be severely 
reduced at the low energies in which this oscillation can occur within a reactor complex. 
These are just the low energies which are necessary for the determination of the hierarchy, 
and so this interference poses a serious problem. One solution is to place the detectors 
orthogonal to arrays of reactors. 

The second interference effect arises from the fact that while the 1-2 oscillation maximum 
is the most useful energy range at which to determine the hierarchy, it enjoys a much lower 
neutrino flux than the 1-2 oscillation minimum. As a result at these energy ranges one can 
expect serious contamination from distant reactors at their 1-2 minimum. This problem, as 
well as the degeneracy between the neutrino flux from distant reactors and ^12, is optimally 
solved by using two detectors at different baselines. However a cheaper solution to the same 
problem is to use a single detector at a shorter baseline, such as 45 km NNW of Daya Bay 
under BaiYiinZhang or, 10 km further west, under the mountains 52 km from Daya Bay. 
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